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Abstract 

We investigate the prospects of detecting a Higgs pseudoscalar (A ) in association with a Z 
gauge boson produced from bottom quark fusion (bb — » ZA°) at the CERN Large Hadron Collider 
(LHC). A general two Higgs doublet model and the minimal supersymmetric standard model are 
adopted to study the discovery potential of pp — > ZA° — > llbb + X, t = e or /x, via bb —* ZA° 
with physics backgrounds and realistic cuts. Promising results are found for m\ < 400 GeV in a 
general two Higgs doublet model when the heavier Higgs scalar (H°) can decay into a Z boson and 
a Higgs pseudoscalar (H° — > ZA°). We compare the production rates from bottom quark fusion 
(bb — > ZA°) and gluon fusion (gg — > ZA°), and find that they are complementary processes to 
produce ZA° in hadron collisions. While gluon fusion is the major source for producing a Higgs 
pseudoscalar associated with a Z boson at the LHC for tan f3 < 10, bottom quark fusion can make 
dominant contributions for tan (3 ^ 10. 
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I. INTRODUCTION 



In the Standard Model (SM), the Higgs mechanism requires only one Higgs doublet to 
generate masses for fermions and gauge bosons. It leads to the appearance of a neutral CP- 
even Higgs scalar after spontaneous electroweak symmetry breaking. The LEP2 experiments 
have established a lower bound of 114.4 GeV [l| for the SM Higgs boson mass at 95% 
confidence level. 

A two Higgs doublet model (2HDM) @ has Higgs doublets 0i and 02 with the vacuum 
expectation values v i and V2- There are five physical Higgs bosons: a pair of singly charged 
Higgs bosons H ± , two neutral CP-even scalars H° (heavier) and h° (lighter), and a neutral 
CP-odd pseudoscalar A . The couplings of the Higgs bosons to fermions and gauge bosons 
depend on the ratio of vacuum expectation values (tan/5 = t^/^i) and a mixing angle (an) 
between the weak and mass eigenstates of the neutral scalars. 

The supersymmetry between a boson and a fermion preserves the elementary nature 
of Higgs bosons. The minimal supersymmetric standard model (MSSM) [^] requires two 
Higgs doublets 0i and 02 coupling to fermions with weak isospin t 3 = —1/2 and t 3 = 
+ 1/2 respectively, to generate masses for fermions and gauge bosons and to cancel triangle 
anomalies associated with the fermionic partners of the Higgs bosons. The Higgs potential 
is constrained by supersymmetry such that all tree-level Higgs boson masses and couplings 
are determined by just two independent parameters, commonly chosen to be the mass of the 
CP-odd pseudoscalar (m^) and tan (3. The mixing angle an between the neutral scalars is 
often chosen to be negative (— vr/2 < an < 0). The LEP2 collaborations have set a lower 
bound of 91 GeV and 91.9 GeV Q for the rrih and the m^, respectively. 

Extensive studies have been made for the detection of a heavier MSSM Higgs boson 
(0° = H° or A ) at the CERN LHC 0, 0, & S El El El El El El E3 • For tan/3 < 5, 
A — > 77, H° — > ZZ* — > 4£, and A°,H° — > ti are possible discovery channels. In the SM, 
H — > ZZ* — > 4£ offers great promise for 160 GeV < rriH ^ 1 TeV 3 3 3] with physics 
background from qq — > ZZ and gg — > ZZ 0, Il9j |. However, the heavy scalar in the 
MSSM is observable via H° -> ZZ* only for tan/5 < 10 and m A ^ 350 GeV Q. The 
detection modes A -> Zh° -> l + l~Tf [10] or l + l~bb |3,ElEi and H° -> h°h° -> bb^ Q 



may be promising channels for simultaneous discovery of two Higgs bosons in the MSSM. 
For large values of tan/3, 0° -»• fip [lj El El El Q, and 0° -> rf @, El El, Q , are 
promising discovery channels for the A and the H°. In some regions of parameter space, 
the rates for Higgs boson decays to neutralinos (H°, A —>■ X2X2) are dominant and they 
might open up new promising modes for Higgs detection 

In two Higgs doublet models, there are two complementary channels to search for a Higgs 
scalar and a Higgs pseudoscalar simultaneously: (i) A — ► Zh° [3, El E3 with a coupling 
proportional to cos(/3 — a#) and (ii) H° — > ZA° with a coupling proportional to sin(/3 — 
At the LHC, gluon fusion can be a significant source to produce a Higgs pseudoscalar (A ) 
and a Z boson (gg — > ZA°) via triangle and box diagrams with the third generation quarks 
[20L 21 1 . The top quark loop diagrams make dominant contribution to gg — ► ZA° for 



tan/3 < 10 |2C 

In the MSSM and a 2HDM with Model II Yukawa interactions [22j for the Higgs bosons 
and fermions, one Higgs doublet (0i) couples to down-type quarks and charged leptons 
while another doublet (02) couples to up-type quarks and neutrinos. The A°bb coupling is 
proportional to tan (3 and this coupling can be greatly enhanced by a large value of tan (3. 
In addition, the Higgs pseudoscalar does not couple to gauge boson pairs at the tree level. 
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Therefore, A — > bb is the dominant decay channel for tan f3 ^ 10 or for itla ~ mz+rrih ~ 210 
GeV with tan/? ~ 2. Recently, we demonstrated that gluon fusion could be a promising 
production mechanism to detect pp — > — > ££66 + X via — > ZA° for tan /3 ~ 2 and 
m A < 260 GeV 

In this article, the prospects of the search for a Higgs pseudoscalar (A ) associated with 
a Z boson produced are investigated. We study the discovery potential of pp — > ZA° — > 
iibb + X) via bottom quark fusion (bb — > ZA°) at the LHC The production cross sections 
of ZA° at the LHC in a two Higgs doublet model and the MSSM are discussed in Section II. 
The dominant physics backgrounds from production of iibb and W + W~bb are presented in 
Section III. The discovery potential of ZA° — > iibb is discussed in Section IV. Conclusions 
are drawn in Section V. 



II. THE PRODUCTION CROSS SECTIONS 



We calculate the cross section for pp —>■ ZA° + X via bb — > ZA° in a two Higgs doublet 
model and the minimal supersymmetric standard model with Model II Yukawa interactions 
for the Higgs bosons and fermions. The parton distribution functions of CTEQ6L1 24] 



are employed to evaluate the cross section for pp —>■ ZA° — > iibb + X with the Higgs 
production cross section a(pp —>■ ZA° + X) multiplied by the branching fractions of Z — > ii 
and A —>■ bb. In the Yukawa couplings of <p bb (</>° = A , H°, h°), the bottom quark mass is 
chosen to be the next-to-leading order (NLO) running mass m^m^) J25[, which is calculated 
with mb(pole) = 4.7 GeV and the NLO evolution of the strong coupling j2(j|. 

In Figure 1, we present the Feynman diagrams for bb — > ZA°. The s-channel diagrams 
contain the heavy Higgs scalar H° \qazh oc sin(/3 — an)} and the light Higgs scalar h° 
[gAZh oc cos(/3 — an)} in the intermediate state. Therefore, this discovery channel might 
provide a good opportunity to measure the couplings of ZH°A° and Zh°A°. The t and u 
channel diagrams are proportional to tan/3. When the t and u channel diagrams dominate, 
the cross section of bb — > ZA° can be enhanced by a large value of tan 2 f3. We have checked 
the relevant couplings with the unitarity condition for the total amplitude. At very high 
energy, the longitudinal polarization vector (e^) of the Z boson can be expressed as the 
momentum vector over Z mass (p^/Mz)- Then unitarity requires cancellation among the s, 
t, and u channel diagrams at high energy. 

Figure 2 presents the cross section of pp — > ZA° — > iibb + X via bb — > ZA° as a function 
of tan f3 in a general two Higgs doublet model and the minimal supersymmetric standard 
model. We have chosen m H = rriA + 100 GeV, rrih = 120 GeV, and an = — 7r/4 for the 
general 2HDM. It is clear that the cross section in a 2HDM can be significantly larger than 
that in the the MSSM when the H° can decay into ZA° with m# > + Mz since m# 
and an are free parameters in a 2HDM. In the MSSM with tan/3 > 10, rriA and rrih are 
very close to each other for rriA ~ 125 GeV, while niA and run are almost degenerate when 



tua ~ 125 GeV ll|. Therefore, the decay H° — > ZA° is kinematically inaccessible in the 
MSSM. The cross section for gluon fusion alone is also presented in this figure. While gluon 
fusion is the major source for producing ZA° with tan/3 < 10, bottom quark fusion can 
make dominant contribution for tan (3 ^ 10. 

The Higgs production rate via bottom quark fusion is very sensitive to the choice of 



factorization scale (fip) III, H, 2^. In Table I we present the cross section of pp — > ZA° 



[° — > 

iibb + X via bb — > ZA° in a two Higgs doublet model with tan (3 = 10 and several values 
of rriA for fi F = M z + rriA, {Mz + m A)/2, (Mz + ttia)/^. In our analysis, we have chosen 
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fip = Mz + rriA with the running bottom quark mass mb{mA) in the Yukawa couplings. Our 
numerical values of cross section are comparable to the recent next-to-leading order QCD 
predictions for Z°A° associated production |30j]. 

TABLE I: The cross section in fb without cuts for pp — > ZA° + X —> iibb + X via bottom 
quark fusion (bb — ► ZA°) at = 14 TeV. We choose three values of the factorization scale (hf) 
for tan/3 = 10 and several values of tjia in a two Higgs doublet model with = 120 GeV, 
niH = m-A + 100 GeV and an = — 7r/4. 



H F \mA (GeV) 


100 


200 


300 


400 


500 


M z + m A 


28.3 


1.81 


8.72 x 10~ 2 


9.19 x 10~ 3 


1.92 x 10~ 3 


(M z + m A )/2 


22.3 


1.52 


7.58 x 10- 2 


8.19 x 10~ 3 


1.74 x 10~ 3 


(M z + m A )/4 


16.1 


1.19 


6.25 x 10" 2 


6.99 x 10" 3 


1.53 x 10~ 3 



To study the effects of the Higgs scalar mixing angle {olh) in a two Higgs doublet model, 
we show the cross section of pp — > ZA° + X — > iibb + X versus an in Figure 3 for tan (3 = 
2, 10, and 50 as well as (a) ttia = 150 GeV and (b) rriA = 400 GeV. Also shown are the 
cross sections in the MSSM for tan/3 = 2, 10, and 50. In the MSSM, an becomes almost 
zero for tan/3 ~ 50. For an < 0, the cross section in a 2HDM is significantly larger than 
that in the MSSM. We include contributions from both bottom quark fusion (bb — > ZA°) 
and gluon fusion (gg — > ZA°). 

For rriA > rriz +Tnh and rriA > rn t +mw, the branching fraction of A — > bb is suppressed 
by A — ► Zh° and A — > ti* with real and virtual top quarks. Therefore, the cross section 
of pp —>■ ZA° — > iibb + X for vtia = 400 GeV is much smaller than that for rriA = 150 GeV. 



III. THE PHYSICS BACKGROUND 

The dominant physics backgrounds to the final state of ZA° — > iibb come from gg — > iibb, 
qq — > iibb, i = e or /x and pp — > W + W~bb + X. In our analysis, we actually evaluated gg — ► 
bi^vbi'v and qq — > bi + vbi~v with dominant contribution frompp — > 1 1 — ► bW + bW~ +X . In 
addition, we also consider backgrounds from — > iigb + X, pp — > iigb + X, pp —>■ iigq + X, 
pp — > + Y, and — >■ j + Y, where q = u,d, s, or c and j = g,q or g. The programs 
MAD GRAPH |3l| and HELAS |32j are employed to evaluate the cross sections for all physics 
backgrounds. 

For an integrated luminosity (L) of 30 fb -1 , we require two isolated leptons with px{i) > 
15 GeV and 1^(^)1 < 2.5 in each event. All jets are required to have pr{b,j) > 15 GeV and 
1^(6,^)1 < 2.5. The 6-tagging efficiency (e&) is taken to be 60%; the probability that a c-jet 
mistagged as a 6-jet (e c ) is 10%, and the probability that any other jet mistagged as a 6-jet 
(tj) is taken to be 1%. Furthermore, we require the invariant mass of the lepton pairs with 
opposite signs to be within 10 GeV of Mz, that is \M t i — Mz\ < 10 GeV to be the signature 
of a Z boson. 

For a higher integrated luminosity of 300 fb -1 , we require the same acceptance cuts as 
those for L = 30 fb -1 , except pr(i) > 25 GeV and pr(b, j) > 30 GeV. The 6-tagging efficiency 
(e;,) is taken to be 50%, and the probability that a c-jet mistagged as a 6-jet (e c ) is 14%. We 
found that the px cuts on leptons and bottom quarks are effective in removing most of the 
SM background, while most leptons from the Z decays and most bottom quarks from the 
Higgs decays survive the pr cuts j23[. 
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In addition, we require that the missing transverse energy ( ]ft T ) in each event should be 
less than 20 GeV for L = 30 ftr 1 and less than 40 GeV for L = 300 fb _1 . This cut on 
missing Et along with the requirement on the invariant mass of lepton pairs ( | M t i — Mz | < 
10 GeV) effectively reduce the background from pp — > W + W~bb + X which receives the 
major contribution from pp — > ti+ X. Our acceptance cuts and efficiencies of 6-tagging and 
mistagging are similar to those of the ATLAS collaboration [l6| . 



IV. THE DISCOVERY POTENTIAL AT THE LHC 

To study the discovery potential of pp — ► ZA° — ► iibb + X at the LHC, we calculate the 
background from the SM processes of pp — *> Itbb + X in the mass window of vha ± AM bb - 
with AM bl = 22 GeV. 

We consider the Higgs signal to be observable if the Na lower limit on the signal plus 
background is larger than the corresponding upper limit on the background p, |33j , namely, 

L(a s + a B ) - N^L(a s + a B ) > La B + N^La~B (1) 
which corresponds to 



(2) 

Here L is the integrated luminosity, as is the cross section of the Higgs signal, and a B is the 
background cross section within a bin of width ±AMjj centered at tua- In this convention, 
N = 2.5 corresponds to a 5a signal. 

In CP-conserving two Higgs doublet models with rriA 3> Mz, a CP-even neutral Higgs 
boson with Standard- Model-like couplings may be the lightest scalar. In this decoupling limit 
|34T |. sin 2 (/3 — an) — > 1 and cos 2 (/3 — a B ) — > 0. We show the cross section with acceptance 
cuts in Figure 4 for pp -> ZA° -> ilbb + X in a 2HDM with m H = m A + 100 GeV, m h = 120 
GeV, and a B = /3 — n/2 as well as the cross section in the MSSM, for L = 30 ftr 1 and 
L = 300 fb _1 . The curves for the 5a and 3a cross sections for the ZA° signal are also 
presented. We include contributions from both bottom quark fusion and gluon fusion. With 
a luminosity of 30 fb -1 , it is possible to establish a 5a signal of ZA° — > llbb for wla ^ 200 
GeV and tan/3 ~ 2 or tan/5 ~ 50. At a higher luminosity of 300 fb _1 the discovery potential 
of this channel is greatly improved for rriA ~ 280 GeV and tan (3 ~ 2 or tan /3 ~ 50. In the 
MSSM, it is difficult to observe the Higgs signal of ZA° since m B ~ uia- In both models, 
if ttla > 250 GeV and tan/5 ^ 7, the branching fraction of A — > 66 is greatly suppressed 
when the Higgs pseudoscalar decays dominantly into and tt* with one of the top quarks 
being virtual. 

In Tables II and III, we present event rates after acceptance cuts for the Higgs signal 
(Ns) from bb — > ZA° — >■ ££66 and the background (N B ) as well as the ratio of signal to 
background N S /N B and N s / 'y/N B in a two Higgs doublet model with tan/? = 10 and 50, 
an = — 7r/4. 

The discovery contours for — > ZA° — > ££bb+X via 66 — > ZA° at the LHC are presented 
in Figure 5 for an integrated luminosity of (a) L = 30 fb _1 and (b) L = 300 fb _1 . We show 
5a contours in the (a^,tan/5) plane for wla = 150, 250 and 400 GeV in a general two Higgs 
doublet model with = 120 GeV and m H = rriA + 100 GeV. In addition, we present the 
curve for the decoupling limit with j3 — a B = 7r/2. For L = 30 fb _1 , it will be possible to 



a s > 



1 + 2JLa B /N 
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TABLE II: Event rates after acceptance cuts for the Higgs signal (Ns = as x L) from bb — ► ZA 
and the background (Nb = <jb x L) as well as the ratio of signal to background Ns/Nb and 
Ns I V Nb in a two Higgs doublet model with tan/? = 10 and 50, an = — tt/4, and m# = + 100 
GeV for an integrated luminosity of 30 fb -1 . 



tan/3 = 10 




m A (GeV) 


N S 


N B 


Ns/Nb 


N s /VNb~ 


100 


196 


1.01 x 10 4 


0.019 


1.95 


200 


14 


2100 


0.007 


0.30 


300 


1 


577 


0.001 


0.03 


400 


< 1 


193 


< 0.001 


< 0.01 


tan/3 = 50 




m A (GeV) 


N S 


N B 


Ns/Nb 


Ns/VN^ 


100 


773 


1.01 x 10 4 


0.076 


7.7 


200 


138 


2100 


0.066 


3.0 


300 


31 


577 


0.054 


1.3 


400 


9 


193 


0.045 


0.62 



TABLE III: The same as in Table I, except that the integrated luminosity is 300 fb . 



tan/3 = 10 




m A (GeV) 


N S 


N B 


Ns/Nb 


Ns/VN^ 


100 


897 


2.62 x 10 4 


0.034 


5.5 


200 


75 


1.02 x 10 4 


0.007 


0.74 


300 


4 


3380 


0.001 


0.06 


400 


< 1 


1170 


< 0.001 


0.01 


tan = 50 




m A (GeV) 


N S 


N B 


Ns/Nb 


Ns/VN^ 


100 


3310 


2.62 x 10 4 


0.13 


20.5 


200 


731 


1.02 x 10 4 


0.07 


7.2 


300 


168 


3380 


0.05 


2.9 


400 


46 


1170 


0.04 


1.4 



discover the ZA° signal for \an\ ^ 0.5 and m A < 250 GeV. The higher luminosity (L = 300 
fb- 1 ) greatly improves the reach for m A up to 400 GeV in a large region of the parameter 
space with \cxh\ ~ 1- 



V. CONCLUSIONS 

Bottom quark fusion and gluon fusion are complementary processes to produce a Higgs 
pseudoscalar (A ) and a Z boson at the LHC. While gluon fusion is the major source of ZA° 
for tan/3 < 10, bottom quark fusion can make dominant contributions to the production of 
ZA° at the LHC for tan/3 > 10. 

We have found promising results for pp — > ZA° — * iibb + X via bb — > ZA° in two Higgs 
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doublet models at the LHC with L = 300 fb _1 for ttia ~ 400 GeV, tan/5 ^ 5, |a.ff| ^ 1, and 
mil = m A + 100 GeV. In the MSSM with itla ~ 125 GeV, tua ~ thh, and the production 
cross section of gg — > ZA° is usually small. 

Gluon fusion (gg — > ZA°) offers great promise for uia ^ 260 GeV and tan/3 ~ 2 |23| . 
The production rate of ZA from gluon fusion at the LHC is suppressed by the destructive 
interference between the triangle and the box diagrams as well as the negative interference 
between the top quark and the bottom quark loops, especially when they are comparable 
with tan (3 ~ 7 [20j . 

In a general two Higgs doublet model, the cross section of bb — > ZA° and gg — > can be 
greatly enhanced when the heavier Higgs scalar (H°) can decay into the Higgs pseudoscalar 
and a Z boson. If we take ran ~ mA in a 2HDM, the Higgs signal will be reduced to the 
level of the MSSM. This discovery channel might provide a good opportunity to discover 
two Higgs bosons simultaneously if the heavier Higgs scalar (H°) can decay into a Z boson 
and a Higgs pseudoscalar (A ). 
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FIG. 1: Feynman diagrams for the signal from bb — > ZAP . 
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FIG. 2: The cross section in fb without cuts for pp -> ZA° + X -> «66 + X at ^/s = 14 TeV, as 
a function of tan /?, for = 150 and 400 GeV, in (a) a two Higgs doublet model with rrih = 120 
GeV, mu = niA + 100 GeV and an = — 7r/4 as well as in (b) the MSSM with rriq = nig = [i = 1 
TeV. We show contributions from bottom quark fusion (bb —* ZA°) and gluon fusion (gg —* Z^4°) 
separately. 
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FIG. 3: The cross section in fb without cuts for pp — ► ZA° + X — > llbb + X at ^/s = 14 TeV, 
as a function of the Higgs scalar mixing angle an, in a two Higgs doublet model with = 120 
GeV, m H = m A + 100 GeV with tan (3=2, 10, and 50, for (a) m A = 150 GeV and (b) m A = 400 
GeV. Also shown are the cross sections in the MSSM for tan/3 = 2 (diamond), 10 (square), and 
50 (circle). We include contributions from bottom quark fusion (66 — > ZA°) and gluon fusion 
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FIG. 4: The cross section in fb for pp — ► ZA° + X — ► Ubb + X versus mA a t \fs = 14 TeV, in a two 
Higgs doublet model with m;, = 120 GeV, mjj = 771,4 + 100 GeV and an = /3 — 7r/2 (the decoupling 
limit), for tan/? = 2 (dashed), 10 (dot-dashed), and 50 (solid). Also shown are the 5a (dashed) and 
3cr (dotted) cross sections for the ZA° signal required for an integrated luminosity (L) of (a) 30 
fb -1 and (b) 300 fb -1 . We have applied the acceptance cuts as well as the tagging and mistagging 
efficiencies described in the text. 
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FIG. 5: The 5a discovery contours at the LHC with an integrated luminosity (L) of (a) 30 fb _1 
and (b) 300 fb _1 in the (aH, tan/3) plane for uia = 150 GeV (medium shading), tua = 250 GeV 
(light shading), and 400 GeV (dark shading) in a two Higgs doublet model with mh = 120 GeV 
and ran = mA + 100 GeV. The discovery region is the part of the parameter space above the 
contours. In addition, we present the curve for the decoupling limit with j3 — an = vr/2. The Higgs 
signal includes contributions from bb — ► ZA° alone. We have applied the acceptance cuts as well 
as the tagging and mistagging efficiencies described in the text. 
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